Charge Ordering and Ferroelectricity in Half-doped Manganites 
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By means of density-functional simulations for half-doped manganites, such as pseudocubic 
Pro.sCao.sMnOs and bilayer PrCa2Mn207, we discuss the occurrence of ferroelectricity and we ex- 
plore its crucial relation to the crystal structure and to peculiar charge/spin/orbital ordering effects. 
In pseudocubic Pro.sCao.sMnOs, ferroelectricity is induced in the Zener polaron type structure, 
where Mn ions are dimerized. In marked contrast, in bilayer PrCa2Mn207, it is the displacements 
of apical oxygens bonded to either Mn^^ or Mn ^ ions that play a key role in the rising of fer- 
roelectricity. Importantly, local dipoles due to apical oxygens are also intimately linked to charge 
and orbital ordering patterns in Mn02 planes, which in turn contribute to polarization. Finally, 
an important outcome of our work consists in proposing Born effective charges as a valid mean 
to quantify charge disproportionation effects, in terms of anisotropy and size of electronic clouds 
around Mn ions. 

PACS numbers: 75.85.+t, 75.47.Lx, 71.15.Mb 
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Introduction. "Improper multiferroics" [1], materi- 
als where ferroelectricity is driven by either spin or- 
dering (SO), charge ordering (CO) or orbital order- 
ing (00), constitute a playground for the physics of 
cross-correlation: the coupling between the different 
orderings and structural distortions is indeed much 
stronger and richer than in conventional covalency-driven 
ferroelectrics.[2| Recently, a complex mechanism of fer- 
roelectricity driven by SO and CO has been proposed 
both in half- doped manganites and in nickelates, sup- 
ported by theoretical studies. [3|-l5| In such systems, it is 
considered that slightly charge-diportionated two mag- 
netic ions form dimers via double exchange interaction, 
resulting in electric dipoles. [6[ In the particular case of 
bilayer-manganite Pr(Sro.iCao.9)2Mn207, it has been ex- 
perimentally suggested that the transition between two 
different CO phases is accompanied by the rotation of 
orbital stripes, in turn related to ferroelectricity. [7] 

This letter is meant to provide insights into cross- 
correlation phenomena and the ferroelectric instability 
in half-doped manganites via first-principles approaches. 
Along these lines, there are two delicate problems casted 
by previous DFT studies: i) it is difficult to unam- 
biguously identify the ground state [j, [sl, |9| between two 
different types of atomic (and related electronic) ar- 
rangements: a) centrosymmetric checkerboard (CB) CO 
pattern of Mn^+ i'^lg^l) ^^^ Mn^+ i'^lg^^g) i^ns and 
b) ferroelectric-active "Zener polaron" (ZP) structure, 
where the two equivalent Mn^-^+ ions are dimerized via 
a "bond-centered" charge; ii) in the CB structure, the 
"local charge" of Mn^+ and Mn^+ ions is in general ill 
defined [10| and a negligible charge separation (in con- 
trast with the nominal valence difference of le) is often 
obtained. Here, we address these two problems by re- 
porting our DFT results on pseudocubic Pro.sCao.sMnOs 
(PCMO) and bilayer PrCasMnsOy. 

Half-doped manganites: phenomenology. In PCMO, as 



in many half-doped manganites, the Mn02 layer shows 
the CE-type AFM configuration, consisting of double 
zigzag FM chains antiferromagnetically coupled both in- 
plane and out-of-plane, which also invokes CO in the form 
of Mn^+ and Mn^+ ions, arranged in a checkerboard pat- 
tern, below Tco=245K (much higher than the magnetic 
transition temperature Tn=175K). Between Tn and Tco, 
there is no long-range magnetic ordering, but likely a 
persisting spin-fluctuating behaviour; we remark, how- 
ever, that the system remains insulating in the presence 
of CO. Non-bonding t2g electrons give rise to a localized 
spin, S* = |, responsible for the AFM super-exchange 
coupling, whereas Cg orbitals, strongly hybridizing with 
O 2p orbitals and consequently producing broad bands, 
are responsible for the double exchange mechanism. The 
half-filled ej subband is a typical example of the coop- 
erative Jahn- Teller (JT) effect, where the JT distortion 
of MnOe octahedrons results in opening of the gap in 
the Cg^ band and in an insulating ground state. The 
crystal structure changes symmetry from Pbnm (at room 
temperature) to lower symmetry in the CO phase, with 
two possibilities: CB-like centrosymmetric Pbnm (alter- 
natively, P2i/m) and ZP-like polar P2inm at low tem- 
perature. The stability of these phases and the ferroelec- 
tric/dielectric properties are discussed later. We took 
CD (charge-disordered) Pbnm structure from Ref. [ll| , 
and ZP-CO P2inm structure from Ref.flJ], consistent 
with previous DFT studies. [9] 

In bilayer-manganite Pr(Sro.iCao.9)2Mn207 as well, 
the MnO bilayer shows the CE-type AFM configura- 
tion, with peculiar different CO transitions. Upon 
cooling, the system was shown to undergo two tran- 
sitions at Tcoi ~370 K and at Tco2 ~310 K (while 
Tn ~153K), where the system shows a transition from 
CD phase to the "COl" phase, forming Mn-e^ orbital 
chains along the 6-axis, and then to the "C02" phase, 
forming orbital chains along the a-axis (which requires 



a 90° rotation of the 00 pattern, cfr FigH]). Interest- 
ingly, the orbital rotation occurred concomitantly with 
a rearrangement of Mn^+/Mn^+ cations, leading to a 
different stacking of the CO-pattern between the bi- 
layers. Although second- harmonic-generation measure- 
ments revealed a non-centrosymmetric state [13|, ferro- 
electric/pyroelectric current measurements could not be 
carried out due to leakage problems. Though truly re- 
markable, the study of Tokunaga et al. leaves many open 
questions (which will be addressed in this paper), such as 
how large the polarization is, if and how the polarization 
is coupled to the SO, CO, and 00, the role of correlation 
effects vs structural effects, etc. 



pseudocubic Pr^ gCa^, 5M11O3 








z~0.35 

b 



!m 




(a) COl: PE 



r-^b \ ^ * ^ * r-^b 

a % Mn3+ • Mn^^t • Mn^+t a 
. Mn4+ C Mn3+I l Mn^+I 

(b) C02: FE 



FIG. 1: CO/00 patterns in (a) COl phase, with CO stripes 
running along a, 00 zigzag lines along b and in (b) C02 
phase, with CO stripes running along 6, 00 zigzag lines along 
a, in bilayer PrCa2Mn207. Inside each bilayer, the inter-layer 
coupling is AFM. The Mn- (apical 0)-Mn angle is denoted 
by 0^^ (0^^) between Mn2+ (Mn^+) ions. The blue arrows 
indicate the direction of electric dipole moment caused by 
the difference between 0^^ and 0^^. 

Computational details. DFT simulations were per- 
formed using the VASP 5.2 code^M] and the PAW 
pseudopotentials [15| within the GGA-\-U formalism |16J 
{U=3 eV and J=0 eV for Mn (i-states[l7[; other values 
of U are also tested), i^-point meshes of (3, 3, 2) and 
(2, 2, 1) were used for the Brillouin zone integration for 
PCMO and PrCa2Mn2 07. For the A-site ordering, we 
consider interlayer Ca-Pr-Ca sandwich stacking along the 
a-direction for PCMO (to do this, the Phnm symmetry 
is kept) and Ca-Pr-Ca sandwich along the c-direction for 
PrCa2Mn207. The Sr doping was neglected, since it is 
expected not to be relevant for our discussion, i.e. the 
one that simultaneously optimizes the energy gain deriv- 
ing from structural/orbital/spin degrees of freedom. We 
succeeded in stabilizing all the CO patterns in the insu- 
lating state, even with [7=0. 

Born effective charges and charge disproportionation. 
In PCMO, when considering the CE-AFM configuration 
and the reduction of symmetry from the CD phase into 
the CB phase, Mn sites split into Mn^+ and Mn^+ sites; 
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FIG. 2: An isosurface and sections in ab plane of spin density 
(blue:up, redidown) of Mn-eJ occupied state at (a) ZP and (b) 
CB in PCMO (zc^O) and at (c) at COl and (d) at C02 phases 
in PrCa2Mn207 (2;:::^0.85). Note the ab frame is shifted and 
rotated for a better comparison. The direction of local polar- 
ization P and of asymmetric Cg electron hopping from Mn"^^ 
to Mn^^ sites are shown by translucent arrows, (e) The same 
isosurface and the section of C02, cut in the be plane. Inset; 
Schematic picture of CB and ZP orbital patterns. Dimerized 
Mn ions are surrounded by dotted ellypses. 

on the other hand, in the ZP phase, the two original Mn 
sites (Mni and Mn2), each split into corner and bridge 
sites along the spin chain (e.g. Mnf and Mnf ), where 
Mni and Mn2 ions are dimerized. When comparing the 
local charge of two different Mn sites, obtained by in- 
tegrating the charge density within atomic spheres, the 
charge separation is smaller than O.Ole, whereas we quan- 



tified difference of the local spin moment between Mn^+ 
and Mn^+ ions as 0.04/iB and 0.33/iB in the ZP and CB 
phases, respectively. In both cases, such tiny differences 
(in contrast to the nominal values of le) are affected by 
large uncertainties, since they depend on the choice of 
the input radii within which the charge/spin density is 
integrated. Due to the strong double-exchange nature of 
Mn-e^ electrons [18|, the Mn-d charge is more like "bond- 
centered" rather than "site-centered" , so that the on-site 
charge based on an ionic picture is ill defined. Luo et 
at. proposed the Mn-d orbital occupancy as an appropri- 
ate indicator of charge separation, giving An=0.17e in 
electron-doped CaMnOs-PJl Actually, the difference be- 
tween Mn^+ and Mn^+ sites in the CB phase is clearly 
seen in the occupied orbital shape, shown in FiglS We go 
one step further by observing that the topology (size and 
anisotropy) of the orbitals should be described by "tenso- 
rial "instead of "scalar" quantities. Therefore, here we in- 
troduce the Born effective charge (BEC), Z*^ = 5Pa/5uf3 
{a^j3 = x^y^z) for each ionic displacement u. The diag- 
onalization of Z* leads to eigenvalues Z^/^^y',^' ^^^ eigen- 
vectors, which correspond to the direction of the response 
of charge density around an ionic site upon displace- 
ment. This represents a practical way to characterize 
the anisotropic nature of pd hybridization along Mn-0 
bonds. For example, the Mn^+ ion in CB phase shows 
/ 3.67 2.94 0.6l\ /7.9l\ 



^Mn3+ 
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where the first (second) component of Z^ is parallel (per- 
pendicular) to the 3z^-r^ orbital lobe in Mn02 planes 
and the third component parallel to the c axis (Figl2] in- 
set). This clearly reflects the bonding nature of Mn Cg 
electrons along the spin chain, consistent with the plot 
of the spin-density (up- minus down-spin components of 
charge density) reported in FiglS) the BEC is enhanced 
along the direction where Mn states shows strong pd hy- 
bridization with neighbouring O p-states. On the other 
hand, the Mn^+ ion in the CB phase shows more isotropic 
Zd components almost parallel to the (a, 6, c) axes (cfr 
TabHj). In the ZP phase, we got anisotropic (i.e. Mn^+- 
type) Z^ both at Mni and Mn2 sites. 



Results: energetics and ferroelectricity. As shown 
in Tab ITU the total energy in Pro.sCao.sMnOa shows 
£^CB>£^zp (with the experimental structure), as con- 
sistent with previous DFT study. [9[ The calculated 
polarization (in Tab. |TT1) shows a small a com- 
ponent, Pa=— 3.3/iC/cm^, and a large b component, 
P6=5.7/iC/cm^, originating from two different origins as 
reported in Ref.[9|. The Pa originates from the dimer- 
ization of Mn ions along the orbital stripes (coexis- 
tence of bond-_and site-centered CO mechanisms, as pro- 
posed in Ref.[3]), consistent with the derivation from 
point charge model; Pj'^^=-2.1/iC/cm^. The dimeriza- 
tion is reflected by the experimental structure: along 



TABLE I: Diagonalized Born effective charge Z"^ (e) at Mn 
sites. Subscripts denote the direction; a' and b' are paral- 
lel to diagonal axis, (a ± b)/2, as bonding and non-bonding 
direction, respectively. 

PCMO 



CB 



ZP 



Mn3+ ( 7.9.. 1.2,, 3.9e ) ^^^ ^ ^-^-^ ^'^^^ ^'^^ 

Mn^+ ( 3.1. 5.U 5.8e ) ^^^ ^ ^'^^^ ^.5,, 4.3e 

^^^^^^^^^= Mnf ( 7.2./ 1.8,/ 4.2e 

7.3./ 3.7,/ 3.1c 



Mn^ 



PrCa2Mn207 



CGI 



C02 



Mn3+ ( 6.74, 0.89^ 3.78c ) ^^'^ ^ ^"^^^ ^'^^y ^■^^- 
Mn^+ ( 5.28. 3.94, 3.45c ) ^^a+ ^ ^.67. 4.93, 3.49c 



Mn; 



4.22. 4.28, 3.90c 



the spin chains, the Mn-O-Mn bond angle is modu- 
lated by large and small angles, 159.6° (intra-dimer Mn) 
and 154.8/157.8° (inter-dimer Mn). Considering that 
the large Mn-O-Mn bond angle enhances the double- 
exchange coupling, we expect the gravity center of charge 
to shift towards the middle of Mn-Mn dimers rather than 
being located on Mn sites, resulting in a local electric 
dipole, suggested as the origin of improper ferroelectric- 
ity in Lao.5Cao.5Mn03[4]. 

On the other hand, the P5 component is due to artifi- 
cially imposing the CE-AFM on top of the P2inm crystal 
structure which leads P5=0 in nonmagnetic configura- 
tion. When forcing the CE-AFM order, the different in- 
teraction between parallel/anti-parallel spin sites causes 
a Heisenberg-exchange-driven polarization of purely elec- 
tronic nature. 



Total energy difference, A^ = Ezv — Ecb and 



TABLE II: 

32 -^coi (meV/Mn) for PCMO and PrCa2Mn207 
different values of U (eV). The polarization P*°* 



A^ = Eco2 

with 

(/iC/cm^) is also reported. 



PCMO 


PrCasMnsOr 


U 


A^;^ 


ptot 


U 


AE^ 


ptot 



3 
5 


-51.7 
-59.3 
-50.9 


(-3.56, 5.01, 0) 
(-3.28, 5.66, 0) 
(-3.18, 6.06, 0) 



3 
5 


-43.4 
-47.2 

-48.7 


(0, 2.18, 0) 
(0, 1.80, 0) 
(0, 1.55, 0) 



Ferroelectricity in bilayer manganites. In 

PrCa2Mn2 07, the calculated energy (Tab. |TT1) shows 
that the C02 phase is energetically lower than COl, 
irrespective of the chosen value of U. The calculated P 
in C02 phase has a non-zero and sizable y component, 
P^^''y=l.SfiC/cm'^. As already pointed out in Ref-IH, 
the FE polarization is related to the stacking pattern of 
the CO bilayers. As shown in Fig. [H the GdFeOa-like 
tilting of MnOe octahedron largely displaces apical 
O (O^P) ions away from the bond center along the b 
direction. One expects the displacement of O^p ion to 
be suppressed by an elastic energy contribution, when 
Mn-e^ electrons form strong Mn-O^^-Mn bonding along 



TABLE III: Mn-O^P-Mn angle 6 (°) in CD, COl, and C02 
phases. 



^33 


^44 


CD 162.80 


162.80 


COl 166.41 


159.22 


C02 162.91 


162.63(A), 159.28(B) 



the c direction. Indeed, the Mn-O^^-Mn bond angles 0^ 
summarized in Tab. IIIH reveal ^^^ < 0^^^ i.e. O^^ ion is 
more displaced between Mn^+ ions, with respect to the 
one between Mn^+ ions. The coupling between the O^p 
displacement and the staggered GdFeOa-tilting pattern 
causes a local dipole inside the bilayer; the latter cancels 
out with the nearby bilayer in the COl phase (resulting 
in antiferroelectricity ) , whereas it gives rise to a net P 
in the C02 phase. 

Furthermore, we observe that the polarization value 
in C02 is close to the P calculated by using PCM 
PP^^(CD)=1.9 /iC/cm^, considering all Mn3-^+ with 
homogeneous valence and Pj^^(C0)=1.8 /iC/cm^, con- 
sidering Mn^+/Mn^+ mixed valences. The almost iden- 
tical value of PP^M(CD) and P^^^{CO) shows that the 
origin of P lies in the apical O ion displacement, and is 
not crucially dependent on the Mn^+/Mn^+ CO pattern 
per se. In addition, this marks a clear difference with 
respect to purely CO-induced ferroelectricity, occurring 
for example in Fe304 and in iron-based fluorides. 19|, |20| 

Further insights can be gained by looking at the 
electronic structure. The 00 both in the COl and 
C02 phases can be clearly seen in the spin-density 
plot reported in Fig|2l showing a similarity between 
COl and ZP, and between C02 and CB phase in 
PCMO, respectively (Figl2j). In analogy with ZP, 
COl shows a dimerization between Mn ions so as to 
induce polarization in the Mn02 plane; however, the 
net induced polarization is canceled out by inter-layer 
antiferroelectric stacking. The difference of C02 with 
respect to CB is represented by two different Mn^+ sites, 
namely Mn^^ and Mug^, and by the anisotropic orbital 
shape. A careful look at the spin density plot reveals the 
anisotropy of the Cg electron hopping. In the C02 phase, 
the charge on the Mn^^ site shows a polar behaviour, 
slightly deviated from x'^-y'^ orbital shape, elongated 
and pointing toward one of the neighboring Mn^+ sites. 
This is consistent with the dimerization of Mn ions, as 
present in the experimental structure, i.e. the larger 
angle of Mn^+-0-Mn^+ than of Mn^+-0-MnB+ indicated 
in Fig. [21 Although the ZP state is realized thanks to 
the double-exchange nature of Mn eg electrons, here 
the dimerization seems to be induced by an external 
factor, i.e. the effective electric field induced by O^p 
displacements. The spin-density plot in the be plane, 
shown in Fig. Efe), reveals the lobes of Mn charge 
pointing along the e direction, forming the inter-layer 



Mn-O-Mn bonding. Remarkably, in both CO phases, 
Mn ions form dimers so as to induce a local P^^ in 
Mn02 planes, that cooperatively adds to the local P^p 
caused by the apical O^p displacement, showing an 
enhancement of the total P. By using the calculated 
BEC tensors, the local contributions to P from each 
layer can be calculated separately and given by the sum- 
mation of the product of the BEC (Z*) by the relevant 
ionic displacement. We therefore obtained Py{CdiO^'^)=- 
0.36/iC/cm2, P^(PrO^P)=P^P=0.65/iC/cm2, and 

Py{MnOf)=P'^^=0.71iJ.C/cm^ per layer, whereas the 
total P5=1.8/iC/cm^. This notably shows that the 
polarization, P^p and P^^, originating from different 
mechanisms and different ions, give quantitatively 
similar contributions. The clear difference in the 
mechanism driving ferroelectricity between PCMO 
and PrCa2Mn207 appear in Z^ (Tbl. [J): all Mn ions 
show a strongly anisotropic Z^ along dimers in the ZP 
phase, whereas isotropic Z^ of Mn^+ is obtained in the 
C02 phase. This implies that in PrCa2Mn207 the Cg 
double exchange mechanism is not relevant but rather 
structural effects play a role in the ferroelecricity. 

Conclusions. We have investigated the ferroelectric 
half-doped manganites and have made a comparison be- 
tween pseudo-cubic PCMO and bilayer PrCa2Mn207. In 
PCMO, although there is an arbitrariness in the crystal 
structure (which was proposed to be either polar or not), 
a polarization is expected in the ZP structure. In bilayer 
PrCa2Mn207, the intercalation of layers gives rise to a 
peculiar situation, where both the Mn dimerization pro- 
cess and the apical oxygen displacements contribute to 
ferroelectric polarization, predicted to be of the order of 
several /iC/cm^. We finally note BEC tensor is a good 
measure to quantify CO and 00 in manganites. 

Acknowledgments. KY thanks Y. Taguchi, Y. Toku- 
naga and D. Okuyama for fruitful discussions and ac- 
knowledges kind hospitality at CNR-SPIN L'Aquila, 
where the manuscript was written. The research lead- 
ing to these results received funding from the Euro- 
pean Research Council under the European Commu- 
nity 7th Framework Programme FP7 (2007-2013)/ERC 
Grant Agreement no. 203523 and from JST-CREST 
"Creation of Innovative Functions of Intelligent Materi- 
als on the Basis of the Element Strategy" . Some figures 
are plotted using the program VESTA[21i]. 



[1] M. Mostovoy and S.W. Cheong Nature Mater 6, 13 

(2007). 
[2] S. Picozzi, and C. Ederer, J. Phys.: Cond. Mat. 21, 

303201 (2009). 
[3] D.V. Efremov, J. van der Brink, D. I. Khomskii, Nat. 

Mater. 3, 853 (2004). 
[4] G. Giovannetti, S. Kumar, J. van den Brink, S. Picozzi, 

Phys, Rev. Lett. 103, 037601 (2009). 



[5] G. Giovannetti, S. Kumar, D. Khomskii, S. Picozzi, J. van 

den Brink, Phys. Rev. Lett. 103, 156401 (2009). 
[6] D. V. Efremov, J. van den Brink, and D. Khomskii, Na- 
ture Matt. 12, 856 (2004). 
[7] Y. Tokunaga et al.. Nature Mater. 5, 937 (2006). 
[8] C. H. Patterson, Phys. Rev. B 72 085125 (2005). 
[9] G. Cohzzi, A. Fihppetti, and V. Fiorentini, Phys. Rev. 
B 82 140101(R) (2010). 
[10] W. Luo et al., Phys. Rev. Lett. 99, 036402 (2007). 
[11] Z. Jirdk et al., J. M. M. M. 53, 153 (1985). 
[12] A. Daoud-Aladine, J. Rodriguez-Carvajal, L. Pinsard- 
Gaudart, M.T. Fernandez- Diaz, A. Revcolevschi, Phys. 
Rev. Lett. 89, 097205 (2002). 
[13] H. Itoh, et al., Appl. Phys. Lett. 96, 032902 (2010). 
[14] G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 



(1996). 
[15] P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 
[16] V.L Anisimov, F. Aryasetiawan and A.L Lichtenstein, J. 

Phys.: Cond. Mat. 9, 767 (1997). 
[17] Weidong Luo, et. al, Phys. Rev. Lett. 99, 036402 (2007). 
[18] V. L Anisimov, L S. Elfimov, and M. A. Korotin, K. 

Terakura, Phys. Rev. B. 55, 15494 (1997). 
[19] K. Yamauchi, T. Fukushima, S. Picozzi, Phys. Rev. B 

79, 212404 (2009). 
[20] K. Yamauchi and S. Picozzi, Phys. Rev. Lett. 105, 

107202 (2010). 
[21] K. Momma and F. Izumi, J. Appl. Crystallogr., 41, 653 

(2008). 



